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We study the traversal times of electromagnetic pulses across dispersive media with negative di- 
electric permittivity (e) and magnetic permeability (p) parameters. First we investigate the trans- 
port of optical pulses through an electrical plasma and a negative refractive index medium (NRM) 
of infinite and semi-infinite extents where no resonant effects come into play. The total delay time of 
the pulse constitutes of the group delay time and the reshaping delay time as analyzed by Peatross 
et al For evanescent waves, even with broadband width, the total delay time is negative for an 
infinite medium whereas it is positive for the semi-infinite case. Evidence of the Hartman effect is 
seen for small propagation distance compared to the free space pulse length. The reshaping delay 
mostly dominates the total delay time in NRM whereas it vanishes when e(ui) = 

Next we present results on the propagation times through a dispersive slab. While both large 
bandwidth and large dissipation have similar effects in smoothening out the resonant features that 
appear due to Fabry-Perot resonances, large dissipation can result in very small or even negative 
traversal times near the resonant frequencies. We investigate the traversal and the Wigner delay 
times for obliquely incident pulses. The coupling of evanescent waves to slab plasmon polariton 
modes results in large traversal times at the resonant conditions. We also find that the group 
velocity mainly contributes to the delay time for pulse propagating across a slab with refractive 
index (n) = -1. The traversal times are positive and subluminal for pulses with sufficiently large 
bandwidths. 

PACS numbers: 

I. INTRODUCTION 

The time for light to traverse through a dispersive medium is interesting and important. A popular measure for the 
delay time of pulses is the Wigner delay time[4| t w — i.e, the frequency derivative of the phase of the output 

wave evaluated at the carrier frequency uj. The Wigner delay time becomes inaccurate for large pulse bandwidths or 
when there is a large deformation in the pulses. 

For pulses and particularly broadband pulses, Peatross et al. [l[ showed that the arrival time of a pulse at a point r 
can be well described by a time average over the component of the Poynting vector S normal to a (detector) surface 
at r as 

u ■ f™^ tS(r, t)dt 

Wr = Too ^7 T\ 77 • (1) 

Here u is taken to be the unit vector along the normal to the given surface. The time of traverse between two points 
(ri,rf) is equal to the difference of the arrival times at the two points, and was shown analytically to consist of two 
parts: a contribution by the spectrally weighted average group delay at the final point Yf 

u • 8(T f ,w) P Rek/dco) • Ar] dco 
A * G = r°° g , rj , {*) 

and a contribution that could be ascribed to the reshaping of the pulse 

At R = T [exp(-lmk ■ Ar)E(r 2 ,u)]-T [E(r l; u)] , (3) 
which is calculated with the spectrum at the initial point r^. Here the operator T is 



J—oo ouj ^ 5 / 



duj 



T\E(t,w)] = 1 r oo " : , (4) 



which represents the arrival time of a pulse at a point r in terms of the spectral fields. The Poynting vector is 
represented by S(r, ui) = Re [E(r, a;) x H*(r, u>)]. Here we take the real parts of the quadratic terms since we use 
complex representation for the fields, i.e., e l ( k r_wt ) for a plane wave. 
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The paper is organised into the following sections: The delay times for evanescent pulses and pulse transport 
through a plasma (infinite and semi-infinite with a boundary) arc presented in Section-2. The delay times for pulse 
propagation in NRM (infinite and semi-infinite with a boundary) are discussed in section-3. The traversal times across 
a dispersive slab are discussed in section-4 and we conclude in Section-5 with a discussion of our results and their 
implications. 



II. ARRIVAL TIMES FOR EVANESCENT WAVES 



We will consider the arrival times for pulses composed entirely of evanescent waves. This is analogous to quantum 
mechanical tunneling of a particle under a barrier. Such situations arise directly in the transport of radiation across 
a metal slab or under conditions of total internal reflection. Now consider the complex wave-vector in a medium, 



k 2 = £^-^r- (5) 
In the limit of small imaginary parts of e and fi, one can write 

k r = Re(fc) ~ \/e r fj, r — £iHi— (6) 

ki = l m (k)c-^p±^^, (7) 

where the subscripts r and i indicate the real and imaginary parts of the quantities. Thus for propagating waves, the 
real part of the wave vector depends primarily on e r & jjb r while the imaginary part is directly proportional to £j & /i, 
or the dissipation. This however becomes different for evanescent waves. To make clear the discussion for evanescent 
waves, we will consider an absorbing electric plasma with e r < 0, £j > and [i = [i r . Now, 

(8) 

ki ~ yj\e r \n r -. (9) 

c 

Thus the real part of the wave- vector depends on the levels of dissipation in the medium (ej) and the imaginary part 
of the wave- vector which determines the decay of the wave depends on \e r \. This implies, in turn, that the definitions 
of the group delay time and the deformation delay time given by Eqns. (2) and (3) respectively get interchanged for 
evanescent waves. This is an important difference for the arrival times of evanescent waves from that of propagating 
waves. 




A. Pulse traversal in an unbounded plasma 

Now we will investigate the arrival times for evanescent pulses inside a plasma. For our calculations we will use a 
dielectric medium with the relative dielectric permittivity of the plasma to be of the causal form 

= 1 - , — . r , (10) 

Lo(uj + lj) 

and the relative magnetic permeability fi = l. 

Now consider a pulse of light initially at r t whose electric field in time is given by 

t 2 

E(n,t) = xE cxp[-— ]exp(-iwi). (II) 

We note that the product Cot, where cZ> is the carrier frequency and r is the pulse duration in Eq. (11) gives a criterion 
for broadband or narrowband pulses. Typically we take Cot = 1000 or 100 for narrowband pulses and Cot — 10 for 
broadband pulses. 

Here we consider that both the source of the radiation and the detector are embedded inside the plasma and that 
the plasma is unbounded. This is to avoid any effects of scattering from the boundaries and study the inherent effects 
of the plasma on the traversal times. The distance between the source and the detector is taken to be Ar. 
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FIG. 1: The total delay time, the reshaping delay time and the group delay time represented by -, A and x as a function of 
the carrier frequency Q (a)Delays for narrowband pulses (pr = 100) and Ar =zc/1007 in an unbounded plasma, (b) Same 
as (a), but for a broadband pulse (pr — 10). Frames (c) and (d) correspond to a semi-infinite plasma with the respective 
parameters similar to (a) and (b). The vertical line drawn at w/7 = 10 separates the propagating waves from the evanescent 
waves. Note that the group delay and the reshaping delay times interchange their roles for the evanescent waves in comparison 
to the propagating waves. The straight line going across the graphs denoted as LL, is the light line for free space propagation 
(At = Ar/c). 



In Figs. 1(a) and 1(b), we plot the delay times obtained for pulse traversal inside the unbounded plasma with a 
plasma frequency uo p = IO7. First of all, we note that the total delay time is negative for a large range of frequencies 
and mostly superluminal below the plasma frequency. The total delay time is also dominated by the reshaping delay 
time at frequency u> less than ui p . Inside an infinite plasma, the absorption determines the energy transport, and hence 
the delay time as well. For larger bandwidths, the contribution of the group delay becomes appreciable at higher 
frequencies (ui ~ 77 to IO7 in Fig. 1(b)). We also note that the qualitative behaviour of the delay time does not 
change appreciably with increase in distance between the source and the detector except for the difference in scales 
1- 



B. Pulse traversal through a bounded plasma 



Here we consider the situation of a plasma with a semi- infinite extent. We consider the source of the radiation is 
inside the plasma at a distance Ar from the planar interface with vacuum. The detector is taken to be in vacuum just 
outside the interface. This is more physical because there would be an interface (impedance mismatched) involved 
with the detector anyway. There is not coupling to any surface plasmon modes of the plasma-vacuum interface. The 
interface has an important effect of changing the amount of energy that reaches the detector via the transmittance of 
the interface given by 

T P M = J'J ^ (12) 



where T p (uj) is the Fresnel transmission coefficient for the P- polarized light. £/ = 1 is the relative dielectric permit- 
tivity in free space. Thus the field at the detector is given by 

H(rj, w) = T»e lk - Ar H(r,, w), (13) 

where k is the wave vector in the plasma. 

In Figs. 1(c) and 1(d), we plot the delay times for pulses from a source within a plasma for different bandwidths 
and with a source to interface distance Ar = ZC/IOO7. First of all, we note that the reshaping delay time is negligible 
compared to the group delay time even for to < uj p . The group delay dominates the total delay time for small Ar( 
zc/107 and ZC/IOO7). This important difference from the case of an unbounded plasma results because the presence 
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FIG. 2: (a) The various delay times plotted for large source-boundary distance (Ar =icj r f) and narrowband pulses (Qt — 100) 
for a bounded semi-infinite plasma. The symbols are similar to Fig. 1. (b) Same as (a), but for broadband pulses (lot = 10). 
The Hartman effect can be observed in graphs (c) and (d) which correspond to the total delay times for various parameters 
shown. 



of the boundary causes a reflected evanescent wave. Now energy transport is primarily determined by the phase 
difference of the incident evanescent wave and the reflected wave, and not only by the dissipation in the medium. 
Thus the group delay time plays the determining role. Secondly it should be noted that the total delay time is almost 
always positive except near the plasma frequency for narrowband pulses and small Ar (zc/107 and ZC/IOO7). One 
notes that negative delay times result for narrowband pulses (ujt = 100) and small Ar even for propagating waves 
(e > 0) near u> = u p . This is a consequence of the k = mode at e = 0. This negativity goes away for larger 
bandwidths. 

At much larger source to interface distances (Ar —zc/j), the deformation of the pulse contributes appreciably 
to the total delay time. We show the delay times in Figs. 2(a) and 2(b) for Ar =zc/j at different bandwidths. 
The behavior of the reshaping delay time (u> < uo p ) tends to that of the behavior in an infinite plasma while the 
group delay time strongly moderates this contribution to the total delay, and the total delay time is positive almost 
everywhere. Surprisingly we note that there is a small region of frequencies where the total delay time goes negative 
even for broadband pulses, and at this large distance involved. However, we note that the spectral width of the region 
where the total time becomes negative, reduces with increasing pulse band width. Hence in the limit of very large 
bandwidths, we expect this spectral width to go to zero asymptotically. 

We note the presence of a Hartman effect in our calculations as well. In Figs. 2(c) and 2(d), we plot the delay 
times with respect to the carrier frequency without scaling with respect to the temporal pulse width. We find that 
over a large range of carrier frequencies below the plasma frequency, the delay time is almost the same for various 
distances involved (Ar =zc/107 and ZC/IOOO7). This is seen for both broadband pulses as well as narrowband pulses 
for these Ar. This saturation of the delay time with distance is a generalization of the Hartman effect for broadband 
pulses which is usually noted for monochromatic evanescent waves. However for much larger distances (Ar =zc/"f) 
which are comparable to the spatial pulse width in free space (I — cr), the deformation takes over and the Hartman 
effect is lost. 



III. ARRIVAL TIMES IN NEGATIVE REFRACTIVE INDEX MEDIUM 



In this section, we will discuss the arrival times for pulses propagating through negative refractive index medium 
with infinite and semi-infinite extent. A sufficient condition for negative refractive index is e < and /1 < at any 
frequency. 

An unusual situation arises when the value of e becomes equal to that of fj, (a case of propagating waves). For this 
case, the reshaping delay turns out to be zero. Note that the value of the refractive index becomes equal to e or \i 
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when e = fi. The electric field at the initial point(in the frequency domain) is given by 

Efc^i^re- 4 ^ 3 (14) 

which is the same field as given by Eq. (11), but calculated in the frequency domain. In an unbounded medium, the 
electric field at the final point is related to that at the initial point as 

E(r / , W ) = E(r i ,u;)e ik - Ar , (15) 

where r/ = n + Ar and k is the wave vector in the medium. The magnetic field is related to the electric field through 
the Maxwell's equation. 

H(r^) =y^— re- 1 ^^, (16) 

and 

Using these, we calculate the delay time for pulse propagation between the initial and the final positions and see that 
the total delay time consists of only one non-zero term which is the group delay time, 

oo r 2 e _ (^ T 2 e _ 2Imk . Ar Rck/duj) . Ar] dbJ 

At = , — . (18) 

roo T2e _(^ T 2 e _ 2Imk . Ar ^ 

J — oo 

This means that the reshaping delay time identically vanishes and no reshaping of the pulse takes place in a medium 
with e — fi. Even for a bounded medium with an interface separating the given medium from vacuum, due to perfect 
impedance matching, the transmission coefficient through the interface is unity. This means that the dispersion in 
the transmission of the pulse through the interface plays no role and the reshaping delay disappears here too. 
We take the Drude Lorentz form for e given by Eq. (10) and a Lorentz dispersion for /x given by 

n=i+- — 4 ■ ■ 

For convenience, we take u>o m = 57, u)^ = QA^ 2 and the rate of dissipation 7 in jj, to be the same as that in e. This 
results in an electric plasma (e < 0, /1 > 0) for < uj < 57, a negative refractive index medium (e < 0, ji < 0) for 
57 < ui < 9.4347, an electric plasma (e < 0, fi > 0) for 9.4347 < uj < IO7 and a positive refractive index medium 
(e > 0, fi > 0) for lo > IO7. 

We study the behaviour in an infinite medium. We plot the delay times(total, group and reshaping delays) for a 
pulse in Fig. 3. We note that at low frequencies, the delay time is negative as in an infinite plasma. However, there 
is a large peak in the total delay time at u> =« 37. The group delay almost exclusively contributes to this. This is 
presumably due to the rapid increase in Re(/i) (as seen from the dispersion curve). In the negative refractive region, 
when there are propagating waves in the medium, there is appreciable contribution from both group and reshaping 
delays. The reshaping delay time is mostly negative. We note that the total delay time is always positive as well as 
subluminal in the negative index frequency region. 

Now, we take a semi-infinite medium with the above material dispersion. As in section 2.2, we take the source to 
be inside the semi-infinite medium at a distance Ar from the interface with vacuum and detect the radiation outside. 
As before, we consider only waves with zero parallel wave vector(k|| = 0). We note that the reshaping delay plays a 
major role in determining the times in the negative refractive index region. In fact for broadband radiation, bulk of 
the total delay time comes from the reshaping delay (Figs. 4(a) and 4(b)). For narrowband pulses (uit = 100), the 
total delay time can become negative near regions where a change of sign occurs in /i or e (Figs. 4(c) and 4(d)) . But 
for broadband pulses (Hit — 10), this negativity disappears even for short distances of propagation. 



IV. TRAVERSAL TIME ACROSS A DISPERSIVE SLAB 



Next we calculate the propagation times for electromagnetic pulses across a dispersive slab, the dielectric permittiv- 
ity (e) and the magnetic permeability (p) of which are given by Eqs. (10) and (19) with respective damping constants 



6 




FIG. 3: The various delay times plotted as a function of Gj in an unbounded medium which behaves like a plasma and with 
negative or positive refractive indices for certain frequency ranges for broadband pulses (lot = 10) and Ar =zc/10y. The 
vertical lines at frequencies 57, 9.4347 and IO7 separate the frequency regions where the field modes are either evanescent or 
propagating. The symbols are similar to Fig. 2. 




FIG. 4: (a) The various delay times plotted as a function of a) in a bounded medium. The parameters and symbols are 
similar to Fig. 3. (b) Same as (a), with Ar =zc/1007. Frames (c) and (d) correspond to narrowband pulses (cut = 100) with 
source-boundary distances Ar =zc/107 and Ar =zc/1007. 



equal to j p and 7„ 



e(u) = 1 



(20) 



MM = — 2—- ' ( 21 ) 

Here uj\ — ljq + ut^ and ^J p ,b,o—^fp,b,o with f p , and fo taken to be 12 GHz, 6 GHz and 4GHz respectively. The 
slab behaves as a positive refractive index medium (PRM) or right handed medium (RHM) (e > 0, fi > 0) when 
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u> > Lu pi a negative refractive index medium (NRM or LHM) (e < 0, \i < 0) within t^o < w < u)f,, and as a barrier 
(e < 0, /i > 0) elsewhere. We take the source of radiation to be placed in vacuum just outside one boundary of the 
slab and the detector just outside the other one. We have taken same medium (vacuum) on either sides of the slab. 
We denote the sides of the source and the detector respectively by regions 1 and 3 and the slab by region 2. So here 
£\ = Hi = £3 = /U3 = 1. £2, and [i2 are respectively given by Eqs. (20), and (21). We take our initial pulse of the 
form given by Eq. (14). The magnetic field is simply obtained using the Maxwell's equations, 

H(rj,w) = y — - re * , (22) 

Here k z \ represents the wave vector in the first medium where the source is present. 

For the P-polarization, the magnetic field at the detector is related to that at the source via the transmission 
coefficient across the slab. The final magnetic and the electric fields at the detector are given by 

H(r ^> ^m^k"^''^' (23) 

and 

E(r,, W ) = A^cVe-^^TH, (24) 
J 2 v / 2^i£3 

where k z3 represents the wave vector in the third medium. 

Here T(w) represents the transmission coefficient across the slab which is given by, 

T H = I _ r >2 e a, 2 A r > ( 25 ) 

where Ar represents the slab thickness and fc Z 2 represents the wave vector inside the dispersive slab, t, t', and r' 
respectively represent the Fresnel coefficients of transmission, and reflection by the slab interfaces and are given by, 

O kzl O k z 2 k z 2 k z 3 

1 £1 1! £2 „/ £2 £3 



kzi _i_ k z2 ' kz2 1 k z 3 ' k z2 1 kz3 ' 

£1 £2 £2 £3 £2 £3 

Here the unprimed, and primed coefficients stand respectively for the coefficients across the first and the second 
boundaries of the slab. For S-polarization, in the expressions of the Fresnel coefficients, the er's are simply replaced 
by /j,'s. 

Also the Fresnel coefficients relate the electric fields across the interface rather than the magnetic fields. Suffixes 
1, 2, and 3 respectively represent the parameters at the source, slab, and the detector sides as described earlier. 
For convenience we later substitute equal material parameters on both (source and detector) sides of the slab. We 
calculate the delay times for different bandwidths (Narrow and broadbands). 

The Wigner delay time was calculated using Eq. (25), and is given by, 

_ d± _ f^tan(fc z2 Ar) +psec 2 (fc z2 Ar)ff Ar 
TbJ duj l+p 2 tan 2 (fc z2 Ar) ' 1 ' 

where 

k z i€2 _)_ k z 26l 

p _ fcz2£l fczl£2 (27) 

for P-polarization and 

k z ip,2 1 k z2 fii 

p = hilEl hl}^ (28) 



for S-polarization. Throughout our calculations (both for normal and oblique incidence), we have checked that the 
Wigner delay time yields the same result as the traversal time for narrowband pulses with the average energy flow 
method. 
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FIG. 5: (Color online) Scaled total delay time of a pulse with different bandwidths (flr) plotted as a function of the carrier 
frequency (f) across a dispersive slab with large thickness (Ar = 10cm). (a) Solid, and dashed lines show delay times respectively 
for a narrowband, and a broadband pulse across a nondissipative slab of RHM. (b) Same as (a), but across a slab of LHM. (c) 
Delay time for a narrowband pulse across a moderately dissipative (y p = O.Ol^p, jm = O.OlWf,) slab (solid line), and a highly 
dissipative (-y p = 0.1u p , j m = 0.1uj b ) slab (dashed line) of RHM. (d) Same as (c), but across a slab of LHM. 



A. Traversal times for normal incidence 



In this case, the parallel component of the wave-vector is zero and the pulse is normally incident on the slab. So 
there is no coupling with the slab plasmon polaritons. Since k x — 0, the Maxwell's equations can be combined to 
give, 

kl = (29) 

This is independent of whether the slab is of a RHM or a LHM. In Fig. 5, we plot the delay times scaled with the 
free space propagation (— ) versus the carrier frequency for both broadband (lot= 10) and narrowband (Cut— 1000) 
pulses. For narrowband pulses, we refer to Q where the Wigner delay times were calculated for a nondissipative 
slab and it was shown that resonant features appear in the delay time behaviors due to presence of the poles of the 
transmission coefficient (Fabry-Perot resonances). We have taken (Ar = 10cm) as large thickness and (Ar = 1cm) 
as small thickness of the slab relative to the wavelength (2.5cm) of the pulse at the electrical plasma frequency (f p ) 
. First we compute the results for the traversal times of broadband and narrowband pulses through a nondissipative 
slab which is achieved by substituting j p = 7,„ = in the expression of e and ji. Figs. 5(a) and 5(b) respectively 
show the traversal times for pulse propagation inside slabs with positive and negative refractive indices. In both the 
figures, it can be observed that, the features due to the slab resonances get smoothened with an increase in the pulse 
bandwidth. So it is expected that for extremely broadband light, these features might completely disappear. Here we 
note that, the results for the narrowband pulses in Figs. 5(a), and 5(b), are exactly the same as those for the Wigner 
delay times given in [5j. Next we study the traversal times for narrowband pulses propagating through dissipative 
slabs of both RHM and LHM (Figs. 5(c), and 5(d)). To include moderate levels of dissipation in the medium, we use 
7 P = O.Olwp and 7 m = O.Olcjf, and for high levels of dissipation, we use 7 P = 0.1lu p and j m = 0-lwb respectively in 
Eqs. (20) and (21). We see that when a small amount of dissipation is introduced in the medium, the time taken for 
transmission through the slab is less than that taken for the nondissipative case. With increased dissipation in the 
slab, one can also clearly observe that the slab resonant features disappear. 

We have also studied the delay times for a highly dissipative slab of LHM for both narrow and broadband pulses 
(Figs. 6(a), and 6(b)). Fig. 6(a) shows the results for a slab with large thickness (10cm) and Fig. 6(b) shows the 
corresponding results for a slab with a small thickness (1cm). It can be clearly seen that the delay time is very small 
near the magnetic resonance frequency. Then it rapidly increases for large frequencies, and after passing through a 
peak, it gradually decreases. For narrowband pulses, the total delay time near the resonance frequency (uiq) even 
becomes largely negative. Even for broadband pulses, with small thickness of the slab, this negativity in the delay 
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FIG. 6: (Color online) (a) Scaled total delay time as a function of the carrier frequency (f) for a narrowband pulse (solid line), 
and a broadband pulse (dashed line) across a highly dissipative (7 P = 0.1cu p , 7™, = O.lc^t,) slab of LHM, and large thickness 
(Ar = 10cm). (b) same as (a), but with a small slab thickness (Ar = 1cm). 



Source of _ 
Radiation 



Detector 



FIG. 7: A layered structure consisting of five slabs used to study the tunneling of pulses. Region 3 is the dispersive slab with 
thickness Ar with £3, and ^3 respectively given by Eqs. (20), and (21). Regions 2, and 4 are air slabs (£2 — ^2 — 1), and with 
thickness d. Regions 1, and 5 are semi-infinite slabs with £1 — £5 = 25, and y,\ = fj,s = 1. The black circle in Region 1 just 
outside the boundary of slab 2 represents the source, and the black screen in Region 5 outside slab 4 represents the detector. 



time appears near loq although to smaller and smaller extent with increasing thickness of the slab. The anomalous 
dispersion of the refractive index of a medium with high amount of dissipation leads to small/negative delay times 
for broadband/narrowband pulses near the resonance frequency. 



B. Traversal times for oblique incidence 



Next, we study the traversal time of a pulse when it is obliquely incident on a slab, i.e., the parallel wave vector 
(k x ) is nonzero. For this case, the Maxwell equations give, 

k » = ]/ ~ k * = -^Vev-ql, (30) 

where k x — q x ~- 

We consider the incident pulses with either S or P polarizations. For convenience, we have scaled all the frequencies 
with respect to the plasma frequency. So here f p = 1, fo = 0.33, and /& = 0.5. The thicknesses of the slab used 
are Ar = -g- (small thickness), for which k p Ar = OAir, and Ar = 2X P (large thickness), for which fc p Ar = 4ir. The 
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thicknesses are relative to the wavelength (X p ) at the plasma frequency. Here we consider pulses for which most of 
the wave vectors are imaginary making the incident wave evanescent. This is achieved by making the second term 
in Eq. (30) under the square root larger than the first term. We calculated both the Wigner delay time and the 
Energy delay time for such pulses most of whose components are evanescent. While the Wigner delay time can be 
easily calculated using the phase shifts, calculating the traversal time via the energy flow for evanescent pulses is a 
non-trivial problem. This is because the energy flow associated with a single evanescent wave in vacuum is zero. Thus 
the Poynting vector for pulse with all spectral components having evanescent wave vectors is zero at the detection 
point in vacuum. Hence it is not possible to compute the energy traversal time for evanescent pulses if the slab is 
embedded in vacuum. One needs to couple the energy in these systems out to measure the pulses. For this, we make 
an arrangement with layered slabs where the evanescent waves are outcoupled to propagating modes in high-index 
media. 

In this arrangement, we have taken two nondispersive slabs of different parameters kept symmetrically on either 
sides of the dispersive slab making a layered slab structure (Fig. 7) . The first and fifth slabs have semi-infinite extent 
with large relative dielectric permittivity (e = 25) and relative magnetic permeability (fi = 1). The second and the 
4th slabs consist of vacuum with e — 1, and \i — 1 and large slab thickness (d) with corresponding k p d equal to (in). 
The 3rd or the middle dispersive slab has e, and \i respectively given by Eqs. (20), and (21) and small slab thickness 
(Ar) with corresponding k p Ar equal to (0.47r). The source is present in the first medium just outside the boundary 
of the second slab and the detector is placed in the fifth medium just outside the boundary of the 4th slab. The value 
of q x in Eq. (30) is chosen in such a manner that the wavevector is real, making the pulse propagating in first and 
fifth slabs, and imaginary making the pulse evanescent in second and 4th slabs. 

First, we plot the Wigner delay time versus both the frequency and the wave vector in a moderately dissipative slab 
for P-polarization. In a rather uniform landscape of delay times, the resonant conditions for the slab surface plasmon 
polaritons (SPPs) stand out in stark contrast where the magnitude of the delay times are comparatively very large. 
Thus the entire dispersion of the SPPs can be traced out (Fig. 8). There are two distinct plasmon modes corresponding 
to the symmetric and antisymmetric modes whose frequencies tend to at large wave vectors. Similarly two modes 
also appear below the magnetic resonant frequency. For highly dissipative slabs also, such plasmon modes are seen 
for evanescent waves, but with large broadening of the dispersion curves (Fig. 8(b)). 

The surface plasmon features are lost when the thickness of the slab is larger than A p . 

Then we studied the traversal times for evanescent pulses having extremely narrow bandwidth (lut = 5000) using 
the energy transport method with our new arrangement of the layered slab structure (Fig. 7). For this arrangement, 
we plotted the delay times for narrowband pulses versus both the frequency and the wavevector, and analysed the 
results for moderate and large dissipative slabs with both P-polarization (Figs. 8(c), and 8(d)) and S-polarization. 
We see that Figs. 8(c), and 8(d) look almost same as Figs. 8(a), and 8(b). Thus, the energy traversal times are also 
significantly affected at the Surface Plasmon Polariton frequencies. It is worth noting that the traversal times are 
large at the resonant conditions. 



C. Propagation through a slab with n = — 1 

Finally, we consider a slab having unit negative refractive index (n = —1) and surrounded by vacuum (n = +1). 
Negative refractive index of unit magnitude can be achieved at a single frequency for a nondissipative slab. The 
properties of such a slab with (n = —1) are very interesting due to the possibility of designing a perfect lens 0, 0] 
By choosing f p = 1, fo = 0.33, /& = 0.5, we get n — —1 at / = ^= (e—-l, fi=-l). With propagation inside the 
medium, the propagation distance increases by a factor of -^g. Here due to perfect impedance matching, no multiple 
reflections take place. Using the expression for group velocity (v g = the group delay along the direction of 

propagation is given by Gd = - ^ s g = (n + ^f^)- For the particular frequency / = -^=, the second term within 

the bracket in the above expression gives a value equal to 4r. We plot both the Wigner delay time, and the group 
delay time versus q x for propagating pulses for the particular frequency mode described above (Fig. 9). From the 
graph, it can be observed that the delay time gradually increases with q x until q x = 1 (where it becomes infinity). 
It can be seen that the graph feature of the group delay time is very similar to the Wigner delay time. Hence it is 
inferred that the group delay mainly contributes to the total delay occuring during the propagation of a pulse inside 
a slab with n = — 1. 
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FIG. 8: (Color online) (a) The Wigner delay times for pulse traversal with P-polarization as a function of scaled carrier 
frequency (f / fp) and scaled parallel wave vector (q x = ^) across a moderately dissipative (7^ = 0.01u p , j m = O.Olu;;,) slab 

with small thickness (Ar = (b) Same as (a), but for a highly dissipative (7^ = O.lujp, -f m = O.loji,) slab. (c)The traversal 
times obtained by the energy transport method across the layered slab structure where the middle dispersive slab is moderately 
dissipative (7^ = 0.01aj p , 7 m = O.OloJt) in nature with small thickness (Ar = -#). (d) same as (c), but for a highly dissipative 
("fp = O.lajp, 7 m = O.lwf,) slab. The dispersion of the slab plasmon polariton modes of the slab stand out clearly and the 
resonant conditions for these modes are characterized by large Energy delay times. 



V. CONCLUSIONS 



In summary, we have shown that the very definitions of the average group delay time and the reshaping delay time 
for the evanescent pulses get interchanged. We have also shown that in an infinitely extended plasma, the delay time 
is primarily determined by the reshaping delay time and is usually negative. In a bounded plasma when the radiation 
is detected outside the boundary (in vacuum), the group delay time dominates and the total delay times are usually 
positive and subluminal for large enough frequency bandwidths associated with the pulses. We also note the Hartman 
effect in the context of energy transport for evanescent pulses when the source to boundary distance is small compared 
to the free space pulse length. In the case of negative refractive index materials, the total delay times are dominated 
by the reshaping delay times. 

For a pulse traversing across a dispersive slab, we have shown that a high amount of dissipation in the slab material, 
along with large pulse bandwidth, smoothen out the resonant features. We have analysed the reason behind the 
occurrence of small/negative delay times near the magnetic resonant frequency which is a consequence of anomalous 
dispersion of the refractive index of the slab medium. We have also shown that the group delay mainly contributes 
to the total delay across a slab with unit negative refractive index, and surrounded by vacuum. 
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FIG. 9: (Color online)The Wigner delay time (solid line), and the Group delay time (cross symbol) plotted versus scaled parallel 
wave vector (q x — ^fp) across a slab with small thickness (Ar = -^), and unit negative refractive index (n = —1) surrounded 
by vacuum. 
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